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The influence of a distribution of relaxation times is studied in annealing experiments. A two- 
relaxation time model is proposed, which permits the calculation of the distribution of relaxation times 
from the crossover data of Spinner and Napolitano as well as that of Ritland. This model also charac- 
terizes the structure of any glass in terms of two parameters. Thus, quenched equilibrium as well as 
nonequilibrium glass were compared with rate cooled glasses with respect to their behavior upon 
further annealing as well as their conductivity at low temperatures with excellent agreement. 

Borosilicate crown glass was found to have a narrow distribution of relaxation times which is 
particular to associated liquid monomers rather than polymers. The results can be explained in 
terms of a topological model for the distribution of relaxation times. The temperature dependence 
of the viscosity is due to a true activation energy rather than a free volume effect. 
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1. Introduction 

This study deals with the analysis of some annealing 
experiments performed on a borosilicate glass by 
Spinner and Napolitano [l]. 1 They established time- 
index of refraction isotherms for their glass in the 
transformation region. It is well known that the index 
of refraction of glass at room temperature is dependent 
on its thermal history, i.e., the type of heat treatment 
that the parent liquid received as it was cooled through 
the transformation region. The index of refraction 
(N) has been used throughout this study because, of 
the various physical properties which reflect this 
change in structure, the refractive index of glass can 
be most accurately measured. It is assumed through- 
out that the mechanical deformations associated with 
change of temperature are purely dilatations. 

The thermal history of the glass in the transformation 
region is shown graphically in figure 1 where the 
indices of refraction (N) have been plotted versus time 
and temperature. The equilibrium indices (line AC) 
for corresponding temperatures are shown on the left 
in figure 1. The glass at temperature 7\, in equilib- 
rium, will have a constant index of refraction A versus 
time. If the temperature is changed instantaneously 
to some lower value T 2 , the index of refraction gradually 
increases to a new equilibrium value C corresponding 
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1 Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. Plot of index of refraction versus temperature and time. 



to the new temperature To. The time dependence 
B—C of this property can be represented as either 

(a) a linear (i.e., nonamplitude dependent) distribu- 
tion of relaxation times or as 

(b) a nonlinear single relaxation (i.e., at any instant 
there is only one relaxation time, but at subse- 
quent intervals the relaxation time changes with 
the distance from equilibrium) [2J. 

Even though any curve can be represented by (a) or 
(b), a systematic series of such curves cannot be 



231 



represented by either (a) or (b) [3] only. Unfortunately 
none of the previous investigators has been able to 
extract either the distribution of relaxation times or the 
nonlinear effects [3]. 

The purpose of this paper is to show how one can 
find a distribution of relaxation times from such 
annealing experiments. 

The time dependence of the physical properties of 
liquids has also been investigated by measuring the 
frequency dependence of these same properties. The 
time and frequency experiments are equivalent and a 
mathematical procedure will be considered which can 
interrelate both types of measurements. 

2. Pertinent Viscoelastic Theory 

When a distribution of relaxation times is present 
in a fluid, the complex modulus of compression, /£*, 
is frequency dependent and can be written [4] 
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where K {) = \im K* 
w^O, 

/< x = limK* 

a> = the angular frequency, 
T = the relaxation times, 
gtr) = the relaxation spectrum. 

Although this is in many ways a convenient form in 
which to represent material properties, it is seldom 
possible to find a unique spectrum to represent a 
measured response function (either a steady-state func- 
tion of frequency like the above or a transient function 
of time like the annealing measurements). This 
results both from the experimental uncertainty of the 
measured function and the fact that a complete re- 
sponse function is seldom measured. As an example 
of the difficulty of finding a unique spectrum, figure 2 
due to Macedo and Litovitz [5] shows a plot of reduced 
modulus versus cor, as well as two different g(r). First, 
g{r) was assumed to be of the form 



g(r) = (bl(rV^)) exp [-/> 2 (ln (r/r')) 2 ] 



(2) 



which is a symmetric spectrum in the logarithm of r 
centered at r' whose width is determined by the value 
of b. 

The second fit was obtained by using two delta func- 
tions for g(r) such that 
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FIGURE 2. Normalized real part of the modulus versus normalized 
frequency for B2O3 at 550 °C from reference 5. 

<p = (t 2 /ti). ,/2 



Depending on the experiment and subsequent data 
reduction, one can obtain the relaxation time spectrum 
at different thermal conditions. In the above example 
where the modulus was calculated from ultrasonic 
velocity measurements, one has a spectrum corre- 
sponding to a stress relaxation modulus at constant 
entropy and strain. In the annealing experiment one 
is measuring the relaxation of a thermally induced 
strain at constant temperature and stress. Even 
though there is a rigorous mathematical connection 
between the constant stress spectrum (retardation) and 
the constant strain spectrum (relaxation) at either con- 
stant temperature or constant entropy, no simple 
relation between the spectra for the two limiting ther- 
mal conditions is available. 

3. Analysis of Annealing Experiments 

3.1. Thermal History of Crossover 

Spinner and Napolitano [1], following some work of 
Ritland [6] and Kovacs [7], transferred a sample from 
the middle of an approach curve having an index (a), 
see figure 1, to a furnace whose temperature (T x ) was 
so adjusted that the equilibrium value of the index (x) 
was equal to the approach index (a). This change in 
temperature was designated as a crossover. The 
changes in index with time at two of these crossovers 
are shown in figure 3. The thermal history of these 
specimens follows: 

(a) All specimens were held at 585 °C until the struc- 
ture was at equilibrium. (Index of refraction was con- 
stant with time.) 

(b) The specimens were reintroduced into the fur- 
nace set at 470 °C for Crossover A and at 500 °C for 
Crossover B. The curves through the open circles 
represent the values of index of the specimens for the 
respective approach curves. 

(c) The change in index with time was followed by 
frequent measurements. When the index had reached 
a value close to the equilibrium value for a tempera- 
ture of 544 °C, the sample was transferred to the 
furnace set at 544 °C and the subsequent change in 
index followed until it had reached its final equilibrium. 
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FIGURE 3. Index versus time at two selected crossovers. 

Equilibrium temperature for both crossovers is 544 °C. Approach temper; 
rossover A is 470 C, and for crossover B, is 500 °C. Data from reference 1. 



Several different starting points for this final crossover 
curve (solid circles) were selected on each of the two 
approach curves, as shown in figure 3. 

3.2. Effect of Distribution of r on Annealing 

Reexamining figure 1, one observes that upon an 
instantaneous change in temperature from T\ to r 2 , a 
thermal strain of N x ( 7\ ) -N«> ( T 2 ) , (B — C),is applied to 
the sample. We will proceed to study the decay or 
relaxation of this strain. 

As noted in the introduction, one cannot represent 
a series of volume-time curves in terms of any given 
sum of exponentials with given relaxation times and 
weighting factors if the total volume change is very 
large. However, we might well expect, if the total 
range of volume covered by the final (crossover) curves 
is small, that a series of such curves could be repre- 
sented by a sum of exponentials without changing the 
constants. We shall find that the curves analyzed here 
can be so represented for any given final equilibrium 
temperature, even though the various crossover curves 
are started from different points on the same approach 
curve or from points on two different approach curves. 
We will term such behavior "linear." 

The simplest such representation which could 
describe the type of behavior shown in figure 3 is two 
exponential terms with equal weighting factors, which 
may be set equal to 1/2. If we were describing the 
index-time relation between two equilibrium conditions 
(starting with the glass in equilibrium at one tempera- 
ture and following the index as a function of time after 
suddenly changing the temperature to a second value 
which is then held constant), we would write: 

N(t) = N 00 +N(l/2 e-'^+l/2 e~' /T *) 

= 1/2 (Mi(t)+M t (t)) (4) 
where Mi(t)=N r > + N i e- tlT i 

where N l =N 2 = N. 



This formulation can be thought of as describing the 
relaxation (change of volume) to a final equilibrium 
value following a change in temperature (analogous to 
an isothermal recovery function) in terms of two 
separate mechanisms, the first of which relaxes ap- 
preciably faster than the second. N\ and A 2 are 
measures of the total relaxation undergone by the two 
separate mechanisms in reaching equilibrium. If the 
starting point is one of equilibrium (Mi(0) = M 2 (0)) the 
two mechanisms will have to relax the same total 
amount, A x (7\)- AUr 2 ), and A x will equal A 2 . How- 
ever, when starting from a nonequilibrium approach 
curve, Ai and A 2 will depend on the particular approach 
curve being followed, as well as on the actual measured 
index at the crossover point, and in general, will not 
be equal. To represent this situation eq (4) can be 
generalized to: 

N(t)= 1/2 (M 1 (t) + M 2 (t)) = N„ 

+ l/2(A 1 e-^ J + A 2 e" r/ ^). (5) 

Figure 4 shows schematically the time dependence 
of A, Mi, and M 2 for a glass initially at equilibrium. 
As the structure rearranges itself to the new equilib- 
rium, Mi will rise faster since it is associated with the 
shorter relaxation time Ti, while M 2 will change more 
slowly. Thus, around the bend of the approach curve, 
there will be the largest difference between Mi and 
M 2 in the glass. As time progresses, Mi will reach its 
equilibrium value for T 2 first and stay there. From 
then on the spread between Mi and M 2 narrows be- 
cause the slow relaxation process continuously ap- 
proaches the same equilibrium value at a slower rate. 
Finally, when the structure is completely at the new 
equilibrium the spread becomes zero, and one has 
again an equilibrium glass (Mi=M 2 ). At a crossover 
(a) the fast relaxation time (r) corresponds to an index 
Mi higher than the (average) measured value, A, and 
the slow relaxation time (r 2 ) corresponds to an index 
M 2 lower than A. 
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FIGURE 4. Illustration of two relaxation times represented by upper 
and lower dashed curves. 

The average or measured index is given by the solid line. The crossover equilibrium 
index x and the transfer point are also shown. 
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FIGURE 5. Crossover A —Solid carve is the average of the two dashed 
carves, each having a single relaxation time; the fast relaxation 
time (ti = 7.77 min). above; and the slow relaxation time (T2 = 64.4 
min) below the eqailibrium line N x = 1.51493. 

Here Vi =0.00103 and A^2 = — 0.00111. Points are third experimental set from cross- 
over A, in figure 4. Data from reference 1. 



If the index at the point of crossover, A(0), is equal 
to the equilibrium index, N(&>) at the crossover tem- 
perature T x , then A 2 = — Ai. A prerequisite for a 
minimum in the crossover is that Ai and A r 2 have oppo- 
site signs. This is fulfilled when N has a value be- 
tween A and B in figure 4. 

The crossover curves were analyzed as follows: 

(a) The logarithm of the difference between the 
equilibrium index and actual index was plotted as a 
function of time. From the best straight line through 
the final points, both r 2 and A 2 were calculated, assum- 
ing t > Ti so M\(t) in eq (5) can be set equal to N x . 

(b) From the initial value of the index at crossover, 
No and the equilibrium value, N x , Ni was calculated 
from 



\ = N n -N x = 



/V,+/V 2 



(6) 



(c) From the time, £ m , at which the index reached 
its minimum value, T\ was calculated from the equation 



dN_ 
dt~ 



:-[At exp (-Wti)]/ti 



-[A 2 exp(-WT2)]/T 2 = 0. (7) 

The results of this computation for the crossover A 
curve with A closest to N x is shown in figure 5. The 



solid line is the average of the upper curve, Mi, and 
the lower curve, M 2 . The agreement of the solid 
calculated curve with the actual experimental points 
is seen to be good. 

3.3. Nonlinearity in Crossover Region 

The definition of linearity adopted requires that the 
relaxation times at a crossover are independent of 
amplitude (i.e., Ni and A 2 ). Since in the linear case 
the relaxation times are characteristic of the equilib- 
rium temperature rather than the approach curve, two 
crossovers having the same equilibrium temperature 
will have the same T\ and r 2 and the same weighting 
factors (V2), regardless of approach temperature. 

The two crossovers in figure 3 present themselves 
as a natural check for linearity. This is done by apply- 
ing T\ and t 2 , already computed for crossover A, to 
crossover B. Equations (5) and (6) were combined 
to give: 



A-A a 



Ni 



[exp (— t/ri) — exp (— */t 2 )] 



+ Aexp(-*/r 2 ). (8) 



Solving for N\ gives: 



M = 2 



(N-N„)-Aexp(-t/T 2 ) 
exp (— t/ri) — exp (— tJT 2 ) 



(9) 



From this equation the average value of N\ for cross- 
over B was obtained using r x and r 2 from crossover A 
and A from eq 6. Figure 6 compares the calculated 
curve (solid line) with the corresponding experimental 
points (the middle set from crossover B in fig. 3). This 
fit is seen to be comparable with that observed for 
crossover A (fig. 5). Linearity of the relaxation process 
in the crossover region with respect to a given equi- 
librium temperature is thus demonstrated to within 
the experimental uncertainty. 

3.4. Further Experiments 

Spinner and Napolitano (1) also reported data for a 
third crossover (C) which corresponds to an initial equi- 
librium temperature of 585 °C, approach temperature 
477 °C, and equilibrium 527 °C. The results of a calcu- 
lation performed as above are shown in figure 7. 
Inasmuch as this crossover contains the largest number 
of experimental points, the agreement of the computed 
curve with these experimental points is even more 
impressive than that for the other two crossovers 
shown. For this crossover, T\ was computed to be 75 
min and r 2 , 540 min. 

An encouraging result is that (t 2 /ti) 52 7°c = 7.2 is 
within the experimental uncertainty of the value of 8.3 
found at 544 °C. Since the distribution of relaxation 
times is temperature independent in most liquids, one 
can assume that this is the case here, even though the 
temperature range is not large enough to justify a 
conclusive statement. 
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FIGURE 6. Crossover B — Solid curve is computed result using fast 
and slow relaxation times found for crossover A. 

In this case /V, = + 0.00069 and ^ = -0.00088. Points arc middle experimental set 
(ref. [ 1 ]) from crossover B, in figure 4. 



1.51590 


1 1 


1 1 


1 1 1 1 1 

EQUILIBRIUM LINE 


1 


1 


1 


80 


1 






c> o 


o 


o 
- — o"~" 


u 

o - 


70 




o/ 








- 


60 












- 


1.51550 


1 1 


1 1 


1 1 1 1 


1 


1 


1 



10 12 14 16 
TIME, HOURS 



20 22 44 26 



Figure 7. Crossover C — Equilibrium temperature 527 °C and 
N x = 1.51584, approach temperature 477 °C. 

Solid line is computed curve based on two single relaxation decay curves (not shown), 
and shows fit to experimental points (ref. [1]). Here T\ = 75 min; N\ = 0.00094: and 
T 2 = 540 min, N 2 = - 0.00086. 



The annealing experiments of Spinner and Napoli- 
tano are based on the following experiment of Ritland 
[6]. A sample of BSC glass was cooled through the 
glass transformation range at a constant rate. Later 
it was reintroduced in a furnace whose temperature 
matched the sample's Active temperature. Figure 8 
shows the time dependence of the specimen after 
reinsertion. Even though the dip is slightly shallower 
the curve is similar. Unfortunately, Ritland did not 
take enough data at "long" times for accurate descrip- 
tion of T2. Since the composition of Ritland's glass is 
very close to that of Spinner and Napolitano, it was 
decided to use the same ratio of T2/T1 = 8 and fit the dip 
by adjusting the average r and the initial spread in 
index. Again, a good fit was obtained, indicating once 
more, that the distribution of relaxation times is a 
property of BSC glass and not of the thermal history. 

3.5. Comparison of t V8 and r s 

The distribution of volume relaxation times has been 
calculated at two temperatures (544 and 527 °C) for 
BSC glass. From these a weighted average relaxation 
time at constant pressure and temperature can be 
defined by: 



T pT = XgiTi = 1/2ti + l/2r 2 



(10) 



On the basis of the values of t p t for each of these tem- 
peratures the apparent activation energy for volume 
relaxation, ///,, was computed from the well-known 
relation 



T pT =Ae H bl RT 



(ID 



where A is an arbitrary constant. Hi, was found to be 
155 kcal/mole (6.49 X 10 5 J/mole). This value is com- 
parable with the activation energy for shear viscosity, 
H s = 140 kcal/mole (5.86 X 10 5 J/mole), over about the 
same temperature range, based on direct measure- 
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FIGURE 8. Solid line is computed curve based on two single relaxa- 
tion decay curves. 

HereT, = 12 min, V, =0.00064 and t 2 = 96 min, /V 2 = - 0.00064, /V x = 1.51453, and shows 
fit to experimental points [6|. Initial distribution of index was obtained by rate cooling 
samples. 



ment of the shear viscosity by the fiber elongation 
method [8]. 

According to a tabulation by Litovitz and Davis [9] 
the ratio of the average adiabatic isochoric relaxation 
time for dilatation, r ls , to Maxwell's average shear 
relaxation time, r s , turns out to be between about 0.79 
and 4.60 for a number of associated liquids. Hence a 
calculation of this ratio should give us a check on 
whether the relaxation times as calculated here to 
represent thermally induced volume changes can 
reasonably be equated to relaxation times associated 
with purely mechanical deformations. 

The isothermal isobaric average relaxation (retarda- 
tion) time, tpt, at 544 °C was calculated from eq (10) to 
be 2.17 X 10 :i sec. It is assumed that the relationship 



tpt(t V 8 = (ktJkt )(C v JC v ) 



(12) 
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derived for a single relaxation mechanism will be ade- 
quate for the comparison desired here. This gives, 
using isothermal compressibilities, kt, and isochoric 
specific heat, C v , from references [11, 12, and 13] 
7^ = 653 sec. 

The average shear relaxation time r s , was computed 
[4] from 

t s = VsIG x (13) 

where G x is the glass modulus of rigidity [14] and at 
544 °C, was found to be 303 sec. Thus, the ratio t vs It s , 
was 2.2, a value in good agreement with the above 
mentioned Litovitz and Davis compilation for asso- 
ciated liquids. Therefore, the relaxation times calcu- 
lated here can reasonably be equated with the relaxation 
process associated with purely mechanical deforma- 
tions. 

The width of the distribution found in BSC glass, 
as given by r 2 /ri = 8 is comparable with that found in 
most liquids. For example, B 2 3 at 550 °C (fig. 2) is a 
little wider, T2/ti = 9, but at higher temperatures it is 
narrower [5]. The associated organic liquids, whose 
distributions are the best known, have about the same 
width provided they are monomeric [9]. For poly- 
meric liquids this ratio is meaningless but would be 
from 10 3 to 10 6 . Thus one can conclude that BSC is 
not polymeric. 

4. Discussion 

The only physical model that would readily predict 
two relaxation times would be two immiscible phases. 
Even though there is a possibility that these relaxation 
times be single they should be separated by many 
orders of magnitude, since one of the immiscible 
phases in the silicate system is usually almost pure 
Si02 with its associated high viscosity (very long re- 
laxation time). Also, had there been such an immisci- 
bility then the temperature cycling of Spinner and 
Napolitano would not have been reproducible. 

We will now present a microscopic model to explain 
the distribution of relaxation times. The topological 
model will consist of a distribution of micro regions 
in space, each with a different relaxation time. Such 
regions have been postulated by various authors. For 
example, Bockris et al. [15], pictured an "iceberg 
structure" for silicate glasses. Litovitz and McDuffie 
[16] and Macedo and Litovitz [5] have proposed that 
the distribution of relaxation times as well as the non- 
Arrhenius temperature dependence of the structural 
relaxation time in associated liquids is due to the co- 
operative behavior of molecules and the possible ex- 
istence of clusters. 

In order to explain the conductivity experiments, 
which will be described below, the glass is assumed to 
have two microdisperse regions which are intercon- 
nected. This should be expected from both the im- 
miscibility work of Haller [17] and the immiscibility 
precursor stage investigated by Ohlberg and Parsons 
[18]. The dimension of these micro regions must be 



smaller than 150 A in order not to be detected by light 
scattering or (ordinarily) by the electron microscope. 
(Zarzycki and Mezard [19] claim to have seen such 
regions even in single component glasses by an elabo- 
rate process using transmission in an electron micro- 
scope.) It is assumed that the relaxation time of 
each region is a function only of its own structure, i.e., 



ri=f(MuT)*f'(Mi 



j*V 



(14) 



This equation is equivalent to assuming that the hy- 
brid equation of Macedo and Litovitz [20] is good for 
each separate region, rather than for the average 
values. 

Evidence for the topological model can be obtained 
from the conductivity experiments of Ritland [6]. 
These conductivity measurements were made well 
below the glass transition and there were no structural 
rearrangements taking place during the measurements. 
The conductivity of a glass below T g is known to be 
due to alkali diffusion through the immobile lattice. 
This ionic diffusion has been shown to have a distri- 
bution of relaxation times. From rate theory one can 
calculate the conductivity due to ionic jumps for a 
single relaxational process to be 



cr- 



ckTri 



(15) 



where a is the "jump distance"; n, the number density 
of ions; T/, the ionic relaxation time; and c, a geo- 
metrical constant. Let us consider the effect of these 
two microdispersed regions on the conductivity. 
First, whatever caused region one to have a shorter 
volume relaxation time, such as higher disorder, less 
intramolecular bonding, etc., may reasonably be as- 
sumed to cause the ionic relaxation time in region one 
to be shorter. Thus, one can write for the conduc- 
tivity of each region: 



<Ti 



ckTru ' 



(16) 



where the nonindex quantities are the same for both 
regions. Since each region is interconnected, the 
measured conductivity can be approximated by 



O" = CT\ -f 0"2 : 



ckT 






(17) 



In fact, if T/, = (l/8)r/ 2 , then the conductivity is con- 
trolled primarily by r/ x . Note that here is where the 
topological model differs from the phenomenological 
approach. Had one assumed a distribution of re- 
laxation times in a homogeneous material then 



(18) 



cA•7 7 (T/ 1 + T/ 2 ) , 



and, provided 8t 7i = r/ 2 , the conductivity would have 
been controlled by r/ 2 . 



236 



Figure 9 due to Ritland shows resistivity versus 
index for glasses that were either cooled at constant 
rate through the glass transformation region (rate 
cooled) or quenched from an equilibrium configura- 
tion (soaked). It can be seen from this figure that 
the variation of resistivity with thermal history for a 
given index is of the same order of magnitude as the 
variation with index for a given type of annealing ex- 
periment. The resistivity is not a single valued func- 
tion of the average structure as represented by the 
index. Thus, any model in which the relaxation is a 
function of the average structure, such as Tool's fictive 
temperature [21], must be ruled out. In the topologi- 
cal model the conductivity is only a function of T/, 
which in turn is a function of M\ and not M 2 . 

In a previous section, it was found that the rate 
cooled glass with an index of 1.51453 had M\ = 1.51517 
associated with T\ and M 2 = 1.51389 associated with 
r 2 . Since as a first approximation resistivity is only 
controlled by structure 1, the resistivity of the rate 
cooled sample having an index 1.51453 should be the 
same as that of the soaked sample having an index 
1.51517 because both have the same values of M\. 
This point is shown plotted in figure 9. If a line is 
drawn through the Ritland's soak specimen data points 
it will intercept his rate cooled curve. This inter- 
section should only occur when the rate sample's 
cooling rate approaches the cooling rate of the chilled 
specimen. Therefore, these lines through the rate 
cooled and chilled specimen data points should be 
almost parallel. Thus for comparison a line is drawn 
parallel to the rate cooled specimen and through the 
"theoretical" point. Since this line is above one 
experimental point and below the other, it is considered 
within the experimental accuracy of the soaked data. 
Note, that eq (18) due to the phenomenological ap- 
proach would predict that o~ is a function of either M 2 
or M\-\-M 2 depending on the assumed functional de- 
pendence of 72, neither of which will fit the data. Even 
though the data are sparse in this experiment it has 




been shown that the two-parameter topological model 
represents the state of the glass regardless of the 
heat treatment. 

Rate Cooling. The best experiment to show whether 
the apparent temperature dependence of the viscosity 
is directly due to temperature or indirectly due to 
structure is the rate cooling process. 

When a glass sample is cooled through the glass 
transition region at a constant rate, initially the sample 
will follow the equilibrium index curve (see fig. 10). 
When the long volume relaxation time, r 2 , becomes 
comparable with the cooling rate, M 2 will lag behind 
the equilibrium value N x . As the temperature con- 
tinues to lower t 2 becomes longer and longer compared 
with the cooling rate, and eventually the structure is 
thermally arrested at a value of the relaxation time, 
tr, which is governed by the cooling rate. Upon 
further cooling M 2 remains constant. M\ will follow a 
similar history and T\ will also be finally arrested at 
the same value, tr. 

Thus as a sample is cooled through the glass tran- 
sition at a constant rate each relaxing process will be 
thermally arrested at the same value of relaxation 
time, tk, which is a function of the cooling rate. If 
the values of the relaxation times are only a function 
of structure, then, at any temperature, all the relaxa- 
tion times in a rate cooled glass cooled at the same 
rate should be equal. The crossover data of Ritland, 
figure 8, clearly demonstrates that this is not the case. 

If we assume that the hybrid equation is completely 
controlled by the activation energy (i.e., /(I*) is given 
by eq (11)), the difference in index M\ — M 2 in the rate 
cooled glass represents the difference in temperatures 
required for T\ to equal r 2 . One, then, can write for 
the activation energy 



//„ = /? In (t 2 /t,){1/7',-1/7 , 2 }- 



(19) 
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FIGURE 9. Resistivity versus refractive index for rate cooled, O, Figure 10. Plot of index versus temperature for a rate cooled 
and soak, 0, samples; X predicted value. Data from reference sample. 

[6\. Equilibrium index (/V^) curve is also shown plotted. 
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in which Ti = T R — N\ 



dT 



and similarly for T 2 ; Tr is 



the crossover temperature of Ritland; and dN/dT was 
taken from Spinner and Napolitano. The activation 
energy calculated from eq (19) is 122 kcal/mol (5.12 
X 10 5 J/mol) which corresponds favorably with the 
140 kcal/mol (5.86 X 10 5 J/mol) obtained from the vis- 
cosity data. 

One must conclude that most of the apparent activa- 
tion energy of viscous flow is due to a true activation 
energy effect. This is in direct contradiction with 
recent papers on the interpretation of the Fulcher 
equation by W.L.F. [22], Cohen and Turnbull [23], 
but agrees with the work of Maurer [24], who an- 
nealed glass under pressure. 

5. Conclusion 

It has been shown how one can obtain the distribu- 
tion of relaxation times by proper analysis of annealing 
experiments. A physical model for the distribution 
was proposed, and it, the topological model, was 
found to explain not only the temperature dependence 
of viscosity, but also the resistivity experiments of 
Ritland. In doing so, Tool's fictive temperature was 
extended to two parameters, which can explain both 
equilibrium and nonequilibrium behavior of glasses. 
The temperature dependence of the viscosity was 
found to be due to a true activation energy rather than 
a free volume effect. The stress (thermal) strain 
experiment gives a distribution many decades narrower 
than the stress (pressure) strain work previously done 
[25, 26] in the relaxation region. Since no pressure 
was applied, no "delayed elastic relaxation" takes 
place in the thermal experiment. Thus, there is the 
added advantage of only measuring the structural 
distribution of relaxation times. 

This analysis was applied to a borosilicate crown 
glass, and it was found that t 2 /ti ~ 8 (Gaussian 6=0.6). 
This narrow distribution indicates that the flow in 
this material is monomeric rather than polymeric 
(10 3 < T2/T1 ^ 10 6 ). The monomeric nature of sili- 
cate glasses has been postulated by Bockris et al. 
[15], from zero frequency compressibility measure- 
ments at high temperatures, but was questioned 
when Westman [27] proved that phosphate glasses 
are polymeric. The present work together with 



Macedo and Litovitz [5] distribution of relaxation times 
in B 2 3 indicates that the iceberg picture of Bockris 
is correct for the borates and silicates. The ratio 
of the average volume to average shear relaxation 
time was found to be ~ 2, in good agreement with 
other associated liquids [9]. 
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